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ABSTRACT 

The Delay Discriminator, an Optimum Tracking Device, is under current 
study as a new space tracking system. The transient analysis of this 
system requires the solution to a non-linear second-order differential 
equation. The Delay-Lock principle is first explained and the implimenta- 
tion of this principle is described. Analysis of the system is discussed 
and a method of digital phase-plane solution is shown. Computer Flow 
Diagrams and Fortran Source Programs are fully described. With out the 
aid of a high speed digital computer, the analysis of this system would 
be difficult, tedious, and inaccurate. 
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1. Introduction. 

The advent of the space age has created many problems concerned with 
position measurements and tracking. In this situation, it is sometimes 
necessary to measure the time delay envolved in sending and receiving a 
given signal when the target is at great distance and moving at high 
velocity. The Digital Delay-Lock Discriminator described in this paper 
is a statistically optimum system which measures the delay between two 
correlated waveforms and satisfies the above distance and velocity re- 
quirements; it mathematically operates in a similar fashion to that of 
the Phase-Lock en The Discriminator is a non-linear system in- 
corporating feedback, and attempts to correlate the delay of the return 
signal with an estimate of the illos generated in the receiver. The 
Delay-Lock is mathematically described by a second-order non-linear 
differential equation, whose solution is not readily obtainable by analog 
methods. The transient analysis of the Delay-Lock Discriminator is the 
main objective of this paper and will be described fully in a later 
section. 

The following section will present a description of the basic 
principles of the Delay-Lock System as first described by Spilker' ^: 
and detail a digital implimentation of this system currently operational 
at the Lockheed Missle and Space Company's Research Labs, This is 
followed by an explanation of phase-plane analysis and methods of graphical- 
ly sketching the transient behavior, 

A further section is devoted to a method of dynamic or sequential 
programming, which allows the computer to step off the operating tra- 
jectories in the phase plane. A CDC 1604 High Speed Digital Computer 


using Fortran Programming provided the transient data which was later 


plotted. 
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This paper then describes a new tracking system applicable to 


many space problems, and presents the transient analysis of its tracking 


operation. 








2. System Operation, 

This section is concerned with the basic theory of operation of 
the Delay-Lock Discriminator as described by Шс. ОЙ A descrip- 
tion of the basic principles of the Delay-Lock System will first be 
presented. After this, a method of generating a pseudo-random sequence 
to represent the signal is described. In conjunction with this signal is 
the techniques of autocorrelation to obtain the desired discriminator 
characteristic. The mathematics of the autocorrelation procedure is 


then completely detailed. All of the above topics are then encompassed 


in a description of the operating Digital Delay-Lock Discriminator. 


2-1 Theory of the DELAY-LOCK DISCRIMINATOR. 

The Block Diagram of the Delay-Lock рінен быс is shown in Fig. 
2-1. The system is a non-linear discriminator incorporating a multi- 
plier, a low pass filter, and a variable delay element. This delay 
element may be ultrasonic, fixed ferrite cores with multiple taps, or 
as in the case that follows, a shift register operating at variable clock 


rates, 
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Fig. 2-1 Block Diagram of the Delay-Lock Discriminator. 
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In an unpublished paper’, Spilker has shown that this discriminator, 

or a slightly modified version is a optimum system in that it provides 

the maximum likelihood (a posteriori, most probable) estimate of delay. 
Defining the delay error ett) such that ELE) = T(t) -T(t) 

where F(t) is an estimate of T(t) and expanding the delayed signal by 


Taylor's Expansion produces the following expression: 


A Р A к 
S(t+T) = S(t ҮТ) +€ S (t rT)+ € s'(t*1)4 Mac 
Ze 
with S(t)normalized to have unity power it can be shown that the output 


of the multiplier X(t) can be given by: 


x ш AB €(+)+nelt) 


where P Belt) is the error correcting term and Ne lt) a noise term. 

In viewing the above equation, the basic principle of the discrimina- 
tor tracking operation is apparent. Consider a sudden increase in the de- 
lay time T(t): Such an E will then increase the initially small 
delay error e (t) . As € (t) increases, the correcting term of the 
multiplier output x(t) increases, and thus the estimate of delay T(t) will 


increase and tend to track the input delay. 


A partially linearized network in Fig. 2-2 describes the multiplier 


FCP) 
| + КАВ F(p)/< 


is linearized when A is constant and describes the locked on region 


output. The close loop transfer function H (P) = 


when € (t)is small. 


LQUIVALENT INPUT 





DEL ^ 
RELAY ч 


KAR Ke) Kinet) ESTIMATE 
LE) KAG Tt) 4+ knelt) | | | 
К КАРТЫ) КАВ 
/1/ 


Fig. 2-2 Partially linearized circuit for delay-lock system. 
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In an unpublished paper’, Spilker has shown that this discriminator, 

or a slightly modified version is a optimum system in that it provides 

the maximum likelihood (a posteriori, most probable) estimate of delay. 
Defining the delay error e(t) such that E(t} = TQ) -T(t) 

where T(t) is an estimate of J(t) and expanding the delayed signal by 


Taylor's Expansion produces the following expression: 


Slt+T)= S(e«T)*€ s'(t rT) +e? s"(t*1)4 -- 2 


With S(t)normalized to have unity power it can be shown that the output 


of the multiplier x(t) can be given by: 


= АВЕ(#)+Пе(+) 


where А Belt) is the error correcting term and Ne lt) a noise term. 

In viewing the above equation, the basic principle of the discrimina- 
tor tracking operation is apparent. Consider a sudden increase in the de- 
lay time T(t): Such an incnegge will then increase the initially small 
delay error e (t) . As e(t) increases, the correcting term of the 
multiplier output x (t) increases, and thus the estimate of delay Та) ап 
increase and tend to track the input delay. 


A partially linearized network in Fig, 2-2 describes the multiplier 
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is linearized when А is constant and describes the locked on region 


output. The close loop transfer function H(P) = 


when é(t)is small. 
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Fig. 2-2 Partially linearized circuit for delay-lock system. 
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Spilker shows in his paper that if s(t)is a stationary (wide 
sense) ergodic random variable with zero mean and slowly varing delays 
A 
те) and 7(%) the loop filter, when optimized, forms the average 


of the multiplier output to obtain: 


E [ke]: E | [As (tT) née) ] «scs f) 
= -KA R; (T-Î) 


Where K. (Q) is the derivative of the auto-correlation function of 
s(t) . Hence М. is not linearly dependent on the delay error but 
functionally dependent upon the error through the differentiated auto- 
correlation function, and thereby causes changes in the effective EN. 
gain. 

By considering the discriminator characteristic for a gaussian low- 
er signal spectrum, it is possible to examine the threshold of interest. 


Such a spectrum and corresponding discriminator characteristic is shown in 


Fig. 2-3, 





Fig. 2-3 Signal power-spectral density and discriminator 
function for Gaussian lowgpass spectrum. ДИ 
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The discriminator characteristic shows only one lock-on region limited by 
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Produces a decrease in delay a ДРЕ here the system is unlocked 
and unstable. Since there is only one region of positive slope, there 
can be no ambiquities in the lock on region and as such, does not possess 
the limitation in tracking that a sine wave signal contains. 

The lock-on performance is more difficult to describe, but when the 
System contains a simple low-pass RC filter, the mathematics are described 
by a first-order non-linear differential equation. When neglecting noise 


terms this equation is as follows: 


--4» 


Je dT,4 = lo) x) --9R, (T-T)/R 


WE Je 


In a later discussion, a simplified form of this equation will be 
used to describe the transient operation of interest. It is well to 
mention here the loop filter characteristics as proposed by Spilker. The 
equivalent filter circuit is shown in Fig. 2-4. The loop filter has 
shown to be optimum for ramp inputs in the presence of white noise, in 
that is minimizes the total squared transient error plus the mean squared 


error caused by noise. 
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Fig. 2-4 Equivalent loop filter and corresponding 
transfer function, /1/ 


The closed loop transfer function for the System is given by: 
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2-2 Generation of the Pseudo-Random Signal. 

This section determines to a large extent the form of the digital 
Delay-Lock System is based on the work of 21 spas Al 

A signal with pseudo-random characteristics and of maximum length can 
be generated with the use of linear feedback shift registers. However, 
not all feedback shift registers will generate maximum length sequences, 
The characteristics of various feedback arrangements are contained in a 
special set of describing polynomials D(x) which are of particular 
concern. 


In considering the binary linear sequential network as shown in 


Fig. 2-5, its operating cycle can be expressed by a T matrix. 


Half-Adder Operation 


1+ 7 = 0 

0 + 0 = 0 

© 14021 

(+) S, Sz у О+1=1 
8 9 =? m — 

S 5, 5; Sy АЖВ + A*B 


Fig. 2-5 Binary linear sequential network with 
four delay elements and half-adder. 


States S So» 5.» 5, describe the internal condition of the network as 


T? 

the four delay element outputs at any time, 54, 5,', S4 » S describe 
the inputs to the delay elements, and will be the outputs of these elements, 
a unit time delay later. Therefore, the Si' are Mod-2 (half adder) sums 


of the S, and the following equations apply: 


SS ES 





Expressed in matrix notation, these equations may be represented as 


follows: 





Or, symbolically, as 
S' = T*S 
where T is called the T-matrix representation of the network. 
Fig. 2-6 shows the state diagram of the present example. The net- 
work generates a maximal length sequence (disregarding the all zero state 


which is trival) representing 2^ - 1 = 15 possible states. 





Fig. 2-6 State diagram of network shown in Fig. 2-5. 
It is not necessary to examine point by point the different states of 
the network to determine if any given arrangement will generate maximum 
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length sequences. This important characteristic can be determined from 
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the network polynomial f(x), where f(x) is defined as follows: 





ф(х) = E - xI 
or in this present example; 


-X 0 1 1 


@(х) = 1 -X 0 0 = х“ + X + 1 
0 1 -х 0 
0 0 1 -х 


In order for a network to generate a maximal length sequence, its 
characteristic polynomial f(x) must be irreducible and not be a divisor 

k ! < n | 
of X- 1 for any interger k Np - 1 where n is the number of delay 
elements. 

The output of any delay element over the complete cycle period M7 sec., 


results in the wave shape shown in Fig. 2-7. 


Ag cu 
Fig. 2-7 Output wave-form from third delay element 
of feedback network shown in Fig. 2-5. 
It is such a sequentíal network that provides the desired signal 
and estimate of delay time in the Digital Delay-Lock Discriminator, 
This network has a characteristic polynomial f(x) as follows: 
ф(х) = х? + x? + 1 
This network is shown in Fig. 2-8, and the output from any one delay 


element constitutes a continuous signal s(t), which has pseudo-random 


characteristics and a maximum length sequence M = 2d - ] or 511 pulses. 








ODD Dp d-D-D р 


D(x) = X94 x 5+ 1 


Fig. 2-8 Maximal length sequential network equivalent to 
linear feedback shift-registers utilized in the 
Digital Delay-Lock Discriminator. 

Although these Sequential Networks have the desirable properties 
for generating a maximal length pseudo-random signal; there is one 
characteristic, the auto-correlation function 14 which has yet to 
be described. The autocorrelation function for the network's sequential 
output signal will be discussed in the following sub-section and shown 


to be ideal, in that it results in a unique and highly desirable Discrim- 


inator Characteristic. 


2-3 Mathematics of the Digital Delay-Lock Discriminator. 

The basic delay-lock discriminator tracking system, covered in 
Section 2-1, can be digitally implemented to track a pseudo-random binary 
signal generated by linear feedback shift register methods previously 
mentioned, The advantage of this new system is that the delay line mention- 
ed in the basic delay-lock system, with its limited tracking range, is re- 
placed by shift registers operating at variable clock rates. 

As mentioned in Section 2-2, a maximum length (2۳ - 1) binary sequence 


generated by the output of any state m of an n stage linear feedback shift 
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register is to be the signal, $. (Е). The amplitude of s (t) will 
range from + 1, with a pulse width ” producing a sequence period of 
м” seconds. Additive noise effects will in general be ignored in the 
following discussion. 

A simplified block diagram of the digital delay-lock network is 
shown in Fig. 2-9. The binary sequence plus noise is fed to a limiter 


network which converts the input signals to binary form. 
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за (£+7) 


Fig. 2-9 Block diagram of the digital delay-lock 
discriminator tracking system. 


The output of the limiter circuit, r(t), is then fed to the primary 

logic circuit which consists of two multipliers (AND gates) and a diode 
resister summing net. The output of the logic circuit, x(t), 15 basically 
tei nary (+1,0,-1) and has an average value which allows the delay-lock 
discriminator to track the received signal. The low pass filter follow- 
ing the logic circuit serves to remove as much of the noise and other 


° 
A 
unwanted signals as possible, The output of the filter, T , is the deriva- 
tive estimate of the delay time, proportional to the radial velocity of 


the target reflecting or returning the transmitted signal. This output 


signal is used to control the clock rate of the receiver shift register. 


11 





Comparison of the clock rates of the transmitter and receiver shift 
registers give a measure of the radial velocity while the Delay is 
measured directly by comparing the instants when both shift registers go 
through a specific state, i.e., the[1,0,0,0,0, anne ‚0,0.|state. Every м7 
seconds a delay measurement is available. 

In general, the theory of operation is very similar to that 
previously mentioned in Section 2-1; that is, the received signal and 
the derivative of the transmitted signal, having been delayed an esti- 
mated time Т, where multiplied to produce an error correcting voltag. 
However, in this case, the received signal is represented by s (t), and 
the delayed version of it represented by s (t m - s. (t -7); both 
signals are then to be multiplied together. The delayed signal is similar 


to the expression of the time derivative of a continuous signal, = in 
that : csi lim s(t+r)- s(&- 7) 
dt ET 

By defining the delay error & (t) to be the difference between the 
True and Estimated delays, e (t)= T e) T(t) , and neglecting 
noise, the output of the multiplier logic circuit, x(t), is defined as 
follows: 

(4,6) = 5 (2.600) [5 66-22 - 5 (6+ 0)] 

Being periodic at intervals MY = 1/ f, seconds, the received signal s(t) 


can be represented by a Fourier series such that: 


Sm (€) =) A sn Sih (nwt) + Åen Cos (nut) 
hco 


Corresponding to this expression, the difference s(t +7) = ES 72) 


can be written as follows: 


vQ 


Smlt+7)-Smlt-?) = -2 ), Б,, 5іһ (nwt) + Beg Cos (nut) 


Heo 
where 


В. = Asn Cos (nu, 7) Ben e ШЕТ Sin (nu. 7) 
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By combining the three above equations, x(t, @ (t)) is as follows: 


xlt et) = Die + уге, есе] 


where D € (t) is defined as the Discriminator Characteristic such that: 


D Lect)] £ P ( Den Sın (n W, €) + Den Cos (nw, €)) 


and nro 
Den = - (Asn Ben- Fkn Ben) ) Den 2 -2 (Asn Bsn + Hen Ben) 

The term ЧЁ, ec) can similarly be expressed, and represents the 
intrinsic noise of the system. 

In considering the Discriminator Characteristic, D( е ) cat 15 
seen that this delay correction term is not explicitly dependent upon t; 
its functional behavior is determined by fixing € and averaging x(t). 
Thus it is apparent that Dle) is the expected value of the multiplier 
output such that: | 
D(e)= E Ir) е) =E Í Sm (t+€)[ Sm(t -7) -Smltr »)] 
This expected value can be expressed in terms of the autocorrelation 
function К ( Q^ ) for the binary sequence s C) so that D( & ) becomes: 

— 

Р (©) = Kmlerr)-Rmle- 7) ; ute ecd) 20 
Thus the Discriminator Characteristic is defined in terms of the auto- 
correlation functions. As was previously stated, the maximum length 
sequence has a unique and desirable R ( С" у; this function is shown 


/5/ 


in Fig. 2-10 . The Discriminator Characteristic is shown in Fig. 2-11: 





Fig. 2-10 Autocorrelation function for maximal 
length binary sequence. 
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Fig. 2-11 Discriminator characteristic for Digital 
Delay-Lock Systems 


It is apparent that D( е ) provides a correcting voltage in a 
linear manner when operating in the region (е|<7? , decreases to zero 


at lel2 27 , and is periodic every M?” seconds. 


2-4 The Digital Delay-Lock Discriminator 
This section describes the operating Digital Delay-Lock Discrimina- 
tor System built under the direction of Dr. James J. Spilker, Jr., at the 
Lockheed Missle and Space Company Research Laboratory, Palo Alto, Calif. 
As shown in Fig. 2-12, the equipment is rack mounted and requires ex- 
ternal analog and digital readout for data analysis. The following major 
components are incorporated in the basic system: 
(a) Rack Mounted Oscillascope 
(b) Transmitter Output Pannel 
(c) Transmitter Shift Register & Logic Circuits 
(d) Receiver Output Pannel 
(e) Receiver Shift Register & Logic Circuits 
(£) Power Control & Voltage Controlled Oscillators 
(g) Loop Filter Inputs and Output, and In-lock 
Tracking Light 
(h) Power Supplies 
The basic logic cards (AND/OR) which comprise the major portion of the 


Transmitter and Receiver sections were scrapped from a defunct Polaris 


Missle check out system and as such do not represent the state of the art 
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in transistor circuitry; however, they provided an inexpensive method 

for realizing the delay-lock system. With the exception of the power 
supplies, the system is completely transistorized; those circuits not 
available originally were built on blank cards. The basic structure of 
the Linear Feedback Shift Registers and associated circuits is shown in 
Fig. 2-14. The primary difference between the two shift registers is 
that the Transmitter shift register operates with a fixed clock rate while 
the Receiver Shift register is controlled by a variable clock rate. The 
clock pulse is generated by a stable crystal oscillator which is fed to a 
doubler or mixer amplifier, This output is fed to a Schmitt Triger cir- 
cuit and then fed to a blocking oscillator. The output is then fed to a 
standard clock-pulse anplifier and Chana the clock inputs of the shift 
register, Although both Voltage Controlled Oscillator networks have a 
manual control for frequency deviation, the receiver VCO is also adjust- 
able by the voltage from the output of the Loop Filter. 

By allowing the difference in clock rates to be large to insure the 
System would not lock-on and track, the error correcting voltage output 
from the filter will conform closely to the system discriminator character- 
istic. When this voltage was fed to an analog plotter, it produced the 
discriminator characteristic shown in Fig. 2-13. This waveform agrees close- 
ly to that predicted for the digital model, the sawtooth discriminator 


characteristic. 


Fig. 2-13 The observed Discriminator Characteristic 
for the Digital Delay-Lock System 
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The outputs of both shift registers were fed to a duel-beam oscilla- 
scope to observe the 1осК-оп апа track conditions; Fig. 2-15 shows 


both binary sequences in an'"unlocked'" and a locked-on and tracking" 


condition. 
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Fig. 2-15 Transmitter and Receiver Binary Sequences; 
(a) in an unlocked condition, 


(b) in a lock-on and tracking condition. 


Fig. 2-16 shows an expanded version of Fig. 2-15 





Fig. 2-16 Expanded version of Fig. 2-15. 
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The block diagram of the complete digital system is shown in 
Fig. 2-17. The system is operable in the presence of large noise 
signals and was tracking smoothly with an input Signal to Noise ratio 
of approximately - 30 db. Spilker has shown the theoretical limit to 


be748db. / !/ 
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3. Phase Plane Analysis. 


This section is concerned with the methods of phase plane analysis 
to obtain the transient behavior of the Delay-Lock Discriminator. The 
first topic to be considered is the basic theory of the phase plane and 
its method of application. Following this will be the development of 
the differential equation which describes the lock-on operation of the 
delay-lock system. The final topic to be considered is the computer 
approach to the phase plane solution for two discriminator characteristics, 
the saw tooth and gaussian cases. 

Appendix A contains the flow diagrams and computer programs written 
in Fortran Source Language for the construction of the phase planes of 


interest. 


3-1. Basic Theory and Application of Phase Plane Analysis. 
A large class of Feedback Control Systems and their operation can be 


described mathematically by a differential equation such as: 


ә 8 

х + POX HK) KH +0 СА 5) = 0 
where P and Q are functions of a signal and its derivative. In phase plene 

| б, | LA 
analysis methods, the value of X is plotted as a function of with the 
5 
individual curves or trajectories for the variation of Я and Ж ап іп- 
dication of time. When examining the transient response of a given system, 
, 
the initial conditions of x and X are plotted and the trajectory (which 
is unique) through this point is the response of the system for these in- 
itial conditions. If the trajectory path converges to a singular point in 
the phase plane, and the system has come to rest with a final value the 
coordinates of the singular point. 
6 
Truxel notes three fundamental restrictions on the use of phase plane 


methods: 
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1, "The phase-plane is useful for the analysis of second-order 
systems Omly......" 

2, "The phase-plane can be useful to study only the transient 
performance of a system subject to initial conditions but other wise un- 
excited." 

3. "The third basic restriction relates to the admissible 
types of non-linearities. The coefficients of x-dot and x can be func- 
tions of x and x-dot, but not of time explicitly." 

Phase plane methods are probably best explained by use of a simple 
example. Consider a second-order differential equation such as: 
о Ф > 
FA tw, XK = ө 


The time operator t can be eliminated by expressing variations in Ж 


e 


| $ 
in terms of A . Define a new variable Y = ж and the equation becomes 
; 2 
+Им ж = O 
Ж 


% , 
Dividing by gives: | dg/ 1 
«X Ч = ہیں ے ال‎ , - ши 2 
№ due X 9 
now define a slope CX such that оќ = 4/2 and again substituting gives: 


X=- Wn № or -Ч/ш = л 
a 


let d Чи or = اا‎ в 4 = = 


The phase plane solution to this equation is found by plotting the iso- 


/6/ 


clines or lines of constant slope. This is shown in Fig. 3-1 along 


with a sample trajectory. 
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Fig. 3-1 Phase plane solution for equation -yl = x/alpha. 
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As shown in the phase plane diagram, the system of concern does 
not converge but circles the origin. Since the equation was that of 
an undamped oscillator, its response should be as shown or simple har- 


monic motion. 


3-2 Application of Phase Plane Analysis to the Problem of Transient 
Response of the Delay-Lock Discriminator. 


The differential equation describing the transient response of the 
Digital Delay-Lock System will now be derived. Figure 3-2 shows a simpli- 


fied version of the system under consideration. 
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Fig. 3-2 Simplified block diagram of the Digital 
Delay-Lock Discriminator. 


The loop filter for this system is a simple low pass RC network as 


shown in Fig. 3-3 along with the filter transfer ration. 





2 
Fig. 3-3 Loop Filter.’ / O Ao > 1 
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The system equat ion can be written in terms of the estimate of the Delay 
A 


derivat ive T , and the Discriminator Characteristic D [ecl such that 
^ 


E or LH | | |+V2'P 
DCE) (p) , wh Ер) is normalized to be Е Ср) = e 
Expressed in terms of operator notation, the equation becomes: 
A 
PT= Fipy> [pce)) 
where p denotes a time derivative operator gt ل‎ 
95” 
ғы the above equation gives: T(t) - V ; Е + is 
Pt: рс) + Рур € 
е до 2 и М 
T * D(9 «Vz D'ce € ся К-Т 
TE = DCE) +Y D'(e)’ € мін 
let X= е Eo 


< i 
-vé « pco«yz О (е): 6 
or -é _ D(é) 
Y +yz DCE) 


For a finite loop gain system, the filter transfer function was 


choosen as FCP) = mL 


system equation becomes: 


P t= Fore): DCE 
[eiF] te Dl + PE De 
THAT =- Dle) V? Dlo E 
уф ка (м6) 027000: = DCE) 
E © _ Dee) -0.1V 
— 7Y-o.tyz De) 


These then are the system equations which when used to construct the 


then 


и 


and in terms of operator notation, the 


phase plane, will describe the transient response of the system. Іп 


order to construct a phase plane that will be of use, it is necessary 
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to evaluate the equation for numerous combinations of slope ү апа 
delay error е . The next topic to be considered describes the use of 


a high speed computer to construct the required phase planes. 


3-3 Construction of Phase-Planes by Digital Methods. 

As previously mentioned, construction of the phase plane for a 
given system can be a time consuming task of point by point evaluation 
to plot the desired isoclines. When the system is a non-linear one, this 
task becomes even more difficult as there may be no obvious symetry to the 
phase plane. Therefore, the problem of evaluating the required isoclines 
to produce transient response trajectories, is simplified when accomplish- 
ed on a high speed computer. Before describing the results of the com- 
puter program, it is first necessary to further define the Discriminator 
Characteristics of interest. 

As mentioned in section 2, the digital system has a saw tooth character- 


istic as shown in Fig. 3-4, and is mathematically defined as follows: 


-6-2 %- -2<64-іІ 


ae - for =| (€ € | 
жи -6 +? tor ILESZ 
О else where 


рсе) Dis, 





Fig. 3-4 Saw tooth Discriminator Characteristic 
The Gaussian Discriminator Characteristic as previously shown in Fig. 
2-3 is also of interest in evaluating the basic delay-lock system. This 


function is mathematically described as follows: 
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The computer programs required to evaluate the isoclines are de- 
tailed in Appendix A; it is only necessary tc say here that many 
hundreds of points were required to complete the phase planes. The 
plot of the isoclines for the saw tooth characteristic is shown in 
Fig. 3-5. It is apparent that the isoclines of magnitude greater than 
5 are concentrated near the е axis and so could not be plotted. A 
sample trajectory is sketched on the phase plane and shows the lock-on re- 
sponse for one set of initial conditions. It is apparent that unless a 
greater number of isoclines is plotted, the accuracy of the trajectory is 
to be deeply suspect. This however is the normal method of solution by 
phase plane techniques. Clearly, this method leaves much to be desired 
and the next section will show two levels of improvement by evaluating 
not only the phase plane but also any specified trajectory, with the 
digital computer. Fig. 3-6 shows the phase plane constructed for the 


Gaussian Characteristic. 
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4. Transient Analysis By Sequential Programming. 


As previously mentioned, the transient response of a given system 
may be obtained by constructing a phase plane associated with the 
system and then sketching the required trajectories by faring the 
curve through the known isoclines. This technique has been recently 
lied to a system known as the Phase-Lock а however, 
this method is at best adequate in that the errors involved in sketch- 
ing are acqumulative, and the final value for lock-on lock critical 
trajectories may be meaningless. Certain problems can be adapted to 
solution by analog computer een, but here again one has to worry 
about amplifier drift during solution producing flucuating errors. For 
the solution of the non-linear equation of the delay-lock, analog methods 
are not practicable, and some other method is desirable. 

This section will present a method of digital analysis by computing 
the desired trajectories with the aid of a high speed general purpose 
digital computer, the CDC 1604. 

Two levels of improvement will be covered; the fixed increment 
approach, and the method of Dynamic or Sequential Programming which 
limits the error to any fixed amount. Actual computer results for the 
transient response of the Delay-Lock Discriminator will be presented, 
Although time is not directly available in the phase plane diagram, it 
is inherently involved and can be obtained by graphical EL 
However, it is available during the computer solution and this method will 
also be covered, 


The Flow Diagram and Fortran Source Language Programs are contained 


in Appendix B. 


4-1 Solution of the Differential Equation by Fixed Increment Method. 


The equation to be solved by phase plane methods as described in 
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Rearranging to solve for the slope y gives: = D(e) Fz D'(e) 
-е 


This is the main value to be computed as the solution of the trajectory 





e 
section 3 is given by: = Е = 


developes. Figure 4-1 shows the basic method of approach. Starting with 
a known Pine in the phase plane representing the initial conditions, B ) E, | 
the slope. By at this point is evaluated. Using a fixed increment of е А 
AE ; the next value of velocity error с is computed in a linear 
€ 
fashion such that é = € tÜ A е and the value of delay at this 
point is determined by с. = Es + AE . This process is continued 
until the trajectory reaches the stable point (in this case the origin) 


where the system is then in lock and tracking, or until the trajectory 


reaches a value of delay error € equal to 2 and the system fails to lock 


in. o 5 e 
€ = =€ +A € 
shpe f Ein 6 
А -LOCK 
| get — N N 
S6 cst To y 
Я AE p — =? 


\42— Lock -ON 





Fig. 4-1 Computer Approach to Evaluation of Trajectory. 
In the method just described, the increment value was chosen to be 
AC = 0.01. Although this method produces satisfactory sets of tra- 
jectories, there is a major difficulty involved, in that when the tra- 


jectory approaches the с axis, the slopes in the phase plane are 
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approaching infinity and even with a very small inciement, the new value 
© 


of е will be made to go too far negative and not be the true trajectory. 
An example of this is shown below and taken from a sample computed trajectory. 


9 
When the value of e first went negative, the routine would set the new 


o 
value of Е 5 @ i » equal to the negative of the old value с апа 


then resume the program. 


Value of £ Value of € Slope 
1.11000 0.29645 -9.07462 
1.12000 0.20569 -13.75419 
1.13000 0.06815 


Switch over occurs as next value 
е ° Е 
of E is negative 


1.13000 -0.06815 1222299 
1.12000 -1.29114 3.83068 


Note this large increase in C in 
just one value of increment, де, 


Figure 4-2 shows the results of this first approach to the evalua- 
tion of the transient response of the system. Several values of initial 


conditions were plotted to describe the lock-on region of interest. 


4-2 Sequential Programming Applied to Phase Plane Analysis. 

Sequential or Dynamic Programming is a basic adaptation of the pro- 
cess developed by Be11man' 19. It is generally described by consider- 
ing the computer as a adaptive device capable of changing system para- 
meters internally as the solution of a given problem develops. In the 
case of the computer solution of the phase plane trajectories, this means 
not operating with a fixed increment, but rather adapting a variable in- 
crement to follow the solution as it develops, 


Instead of the fixed increment as previously covered, consider a 


variable increment AÉ y , Which will take on one value as the slope 
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in the phase plane becomes zero, and decreases to zero as the slope 
approaches infinity. It was decided that the variable increment in 


this problem would be defined as follows: ° 


с = 0.02 
= Bar 


This allows an increment to take on values from 0.02 to zero. When the 





trajectory passes through slope areas approaching infinity, the trajectory 
starts to move almost vertically until this high slope region is passed. 

A comparison of this method with that of the fixed increment type is 

shown in Fig. 4-3. The method of sequential programming then allows the 
trajectory to recover from these high slope regions without being overly 


influenced by them, 





Fig. 4-3 Comparison of fixed and variable 
increment trajectories, 


Figure 4-4 shows the results of variable inciement method for deter- 
mining transient response of the Delay-Lock System, These curves are based 
on a saw-tooth discriminator characteristic and the system having infinite 
loop gain. Various initial velocity error conditions were plotted to 
demonstrate the lock-on characteristics of the system. Considering the 
trajectory describing an initial velocity error of unity, the curve starts 
to increase in velocity error while the system is reducing the delay error 
and correcting the estimate of signal delay. The trajectory starts its 


spiral towards the singular point which corresponds to a lock-on and track 
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condition. In this case the system does lock on and so will tolerate 
a velocity error of this magnitude. However, in the case of an initial 
velocity error of 2,3, the system attempts to lock-in and track but 
does not quite make it and assumes a no-lock condition. The system 
will then remain out of lock until the estimate of the received signal 
is again in the discriminator region. Fig. 4-5 shows the transient 
response for the same system discriminator characteristic but concerns 
the finite loop gain case. The lock-no lock critical trajectory occurs 
for the same initial velocity error but in this latter case, the system 
comes to rest with a finite delay error proportional to the initial 
velocity error. Since a trajectory may be started any where in the 
. phase plane, the computer curves offer a second analysis of the operating 
system. This concerns the case where the system is in lock and tracking 
a target with a constant velocity. If the target accelerates (or de- 
celerates), then the target assumes a new velocity and the system is 
temporarily out of track and must lock on in the manner described above. 

The transient response for the delay-lock system with a Gaussian 
Discriminator Characteristic and infinite and finite loop gains is shown 
in Fig. 4-6 and 4-7. The phase plane curves are expleined in the manner 
previously mentioned, Although the computer data was hand drawn, the 
curves could easily be constructed by a high speed digital plotter to 
give as many trajectories as desired. 

The. error in the trajectory as it passes through the region of large 
valued slopes may be approximated in the following manner: 


e 
The next value of € may be expanded in a Taylor's series to be: 


On 
E (є: Є Се. ғ«Ек)- Со лек ЗЕ + EL AR 


Assuming a parabolic trajectory when crossing mad the infinite slope 


area allows the following simplifications: 
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2-7 
and assuming that е “a М is approximately T Á е 


4-3 Computation of the Lock-On Time Response. 

As previously mentioned, time response is not directly available 
for the phase plane diagram, but may be determined by graphical tech- 
niques. However, this procedure is again a tiresome method which is not 
very accurate. It is possible to compute the time response for a given 
trajectory while the computer routine is in the process of constructing 
the trajectory. Fig. 4-9 shows the typical delay error versus time response 
for several initial velocity conditions along with a typical velocity error 
versus time response. As is apparent, the delay error term r ion. 
at a value of -2 and oscillates in a damped fashion converging to zero 
value delay error when the system has attained Lock-On. These transient 
responses were observed on a analog plotter for the Digital System and 
were of the same general appearance as those predicted. Fig. 4-8 shows 
an observed transient response for the approximate velocity error function 
when the system was in the process of lock-on. In order to show the method 


of computing transient times it is necessary to expand the delay error 


| T los 
function in a Taylor's Series. Therefore, с, (4)- €, AT =. + at е 
(А 


and the increment of time required for the trajectory to pass from 


Es to e is given by: At: EIE дек 
е ® 


> E, 
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It is only necessary to sum the time increments to find the total 


time expired at any given point on the trajectory. The small error 


2 ee 
term ALE was neglected. 


نے 


2! 





Fig. 4-8 Observed wave-form of velocity error 


response of Digital Delay-Lock System 
during lock-on. 
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5. Conclusions. 

The Digital Delay-Lock Discriminator with its pseudo-random signal 
characteristics offers a method of tracking targets at great ranges and 
high velocities, while operating with moderate power requirements. The 
experimental system is now operating at the Lockheed Missle and Space 
Company and is under consideration by NASA as a proposed tracking and 
rondevous system for maned and unmaned satelites. The system has other 
possibilities such as an underwater tracking system. Since the system 
has the capability to track targets when the Signal to Noise ratio is 
іп the order of -30db, the pcwer of the transmitted signal could be 
maintained at a level below that of the ambient noise surrounding the 
target, and as such, would not give the target an indication it was 
under survailence. 

In order that the system be optimized, it was necessary to deter- 
mine the transient response of a normalized system during the lock-on 
process.  Phase-plane analysis methods were then utilized to obtain 
this response. With out the aid of a high speed digital computer, this 
analysis would not have been easily accomplished. The data obtained 
from the sequential programming methods demonstrate that this phase 
plane approach to the analysis problem is applicable to a large number 


of non-linear systems. 
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APPENDIX A 
This Appendix contains the Flow Diagram and Fortran Source 
Programs required to evaluate the isoclines (lines of constant slope) 


of the phase planes which represent the solution to the differential 
equation: с D ( е) 
comm. "n 
Я+үз (©) 





In this problem, it is necessary to fix a value of slope 5 
© 
and evaluate the above equation for Е with Е taking on the range 
-р < EL 2 ; an increment of 0.01 was choosen for : 


Values of slope е used were: + ( 0.2, 0,42 0.6, 0.8, 1.0, 2.055320, 
4.0, 5.0, 10.0, 15.0, 20.0, 25.052300, 2407 40,08 45.0, 50.0 } : 

It was thought that this many values of slope would be ге Е to 
properly bound the phase plane. As it turned out, only slopes in the 
region y < |0 were useful since those of higher value could not be 
plotted on a useable phase plane. The results of these programs are 


shown in Section 3, Figures 3-5 and 3-6, 
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APPENDIX B 
This appendix contains the Flow Diagrams and Fortran Source 
Programs required to generate the desired trajectories describing the 
lock-on behavior of the system. The results of these programs are 


shown in Section 4, Figures 4-2 through 4-7. 
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